Physically Realizable Entanglement by Local Continuous Measurements 
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Quantum systems prepared in pure states evolve into mixtures under environmental action. Phys- 
ically realizable ensembles are the pure state decompositions of those mixtures that can be generated 
in time through continuous measurements of the environment. Here, we define physically realizable 
entanglement as the average entanglement over realizable ensembles. We optimize the measurement 
strategy to maximize and minimize this quantity through local observations on the independent 
environments that cause two qubits to disentangle in time. We then compare it with the entangle- 
ment bounds for the unmonitored system. For some relevant noise sources the maximum realizable 
entanglement coincides with the upper bound, establishing the scheme as an alternative to locally 
protect entanglement. However, for local strategies, the lower bound of the unmonitored system is 
not reached. 
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Decoherence is the process in which the exchange of 
information between a quantum system and an external 
environment continuously downgrades quantum proper- 
ties of the former This dynamics turns initially pure 
entangled states into mixed less entangled ones, destroy- 
ing the capacity of the system for quantum applications. 
In this picture, the description of the properties of the 
system in time results from ignoring the information that 
leaks into the environment and only considering the gen- 
eral statistical characteristics of the reservoir. If, how- 
ever, one is able to recover some of this information by, 
for example, continuously observing changes into the en- 
vironment, the time evolution of the system undergoes a 
different dynamics which is conditioned on the results of 
these observations. Such conditional evolutions are usu- 
ally referred to as quantum trajectories [1] because if the 
system is initially prepared in a pure state, the sequence 
of measurements performed on the environment deter- 
mines a respective sequence of other pure states for the 
system. Each trajectory corresponds to one realization of 
the experiment and the open system evolution is recov- 
ered when averaging over all possible trajectories. Such 
trajectories have already been observed both in massive 
particles and light [2]. 

There are infinitely many different ways to observe a 
physical environment, also referred to as unravellings. 
Each one of them defines a different set of experimen- 
tally realizable trajectories. However, even though at any 
given time the incoherent sum over all possible trajecto- 
ries reproduces the quantum state of the unmonitored 
system, that does not mean that all its pure state de- 
compositions are available, since not all ensembles can be 
achieved through this continuous monitoring process [3]. 
The fact that one can reconstruct the state allows one 
to recover the average value of any observable in time by 
averaging over the trajectories. However, when it comes 
to obtaining the entanglement of the system, the restric- 



tions imposed by the unraveled time evolution may be 
too strong [3]. 

In this paper, we investigate equivalent local mea- 
surement strategies over the independent reservoirs that 
cause two entangled qubits to lose coherence and en- 
tanglement in time. The restriction on the ensembles 
that can be continuously generated in time becomes clear 
when we compare the minimum and maximum aver- 
age entanglement [5] resulting from locally realizable un- 
ravellings to the extremes over all possible decomposi- 
tions of a given density matrix, respectively the entangle- 
ment of formation Ep [6J and the entanglement of assis- 
tance Ea [2] • The locality restriction applied means that 
the measurements performed over independent reservoirs 
should, themselves, be independent as well, or at most 
classically correlated. This local condition is essential 
whenever the decohering qubits are distributed over dif- 
ferent and distinct channels, like in quantum teleporta- 
tion, communication or likewise. 

We analyze both dephasing and spontaneous emission 
as the decoherence sources. For dephasing, we show that 
it is not only possible to continuously reach the entan- 
glement of assistance of the system but the protocol to 
do so does not even require any classical communication 
between the observers, as long as local phase feedback is 
allowed. On the other hand, we also show that the entan- 
glement of formation can never be achieved through local 
observations, even with the addition of classical noise to 
the measurement process [S] , already making explicit the 
impact of the missing decompositions on the realizable 
entanglement. This impact becomes even more explicit 
for dissipative reservoirs where Ea can only be achieved 
for some particular initial states, and, once again, Ep 
can never be achieved through local observations. This 
means that not only some decompositions are not pro- 
duceable through local continuous monitoring of the en- 
vironment, but, more important, neither is the entan- 
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glement of p{t). This also means that the independent 
monitoring of local reservoirs is indeed a good strategy to 
preserve entanglement for quantum communication and 
teleportation protocols, specially in quantum feedback 
schemes 0. 

The physical system we describe is that of two non- 
interacting qubits initially prepared in a pure state a = 
\a) ((t| and coupled to independent reservoirs in the Bohr- 
Markovian approximations [I]. Their unmonitored time 
evolution is determined by a master equation of the type 
p = Lj\,p + LbP where 



/:p^^|(2r,prt-{rtr„p}). 



(1) 



and the set {Tj} represents the subsystem- reservoir cou- 
plings. p{t) can be formally obtained by integrating this 
equation in time. It also corresponds to a completely 
positive map in the Kraus form [lOj acting on the initial 
state. 
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with K representing a vector of Kraus operators K — 
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ith Kraus operator can be identified as a measurement 
operator such that probability is conserved, meaning 
— 1) ^i^d the measurement process is charac- 
terized by a non negative distribution (K^Ki) > 0. The 
K unravelling gives output states representing a decom- 
position of the system state p = ^'5'^^, with K|ct) = 5*, 
and = Ki\a) = jV'i) (the tilde endicates a unnormal- 
ized state). 

This is one of infinitely many possible decompositions 
of p and any other decomposition p = can be ob- 
tained by a unitary rearrangement of — U'^"^ , 
such that p = ^^■f = '^U^U*'^^ = ^f*^ where U is a 
unitary matrix respecting iU^U)^ = U^U* = 1. In [TT], 
the authors show a one to one relation between unravel- 
lings and pure state ensembles by demonstrating that this 
same unitary freedom applies to the unravelling choice, 
or in other words, that the same unitary takes the unrav- 
elling K that generates ^ to another unravelling G that 
generates $, with G^^ — UK^: 



KaK^ = KU^aWK^ = GcrG 



(3) 



Their result shows that in principle all possible ensem- 
bles that decompose the decohered state at a given time t 
can be obtained through measurements on a purification 
of the state. These include, in particular, the decompo- 
sitions that handle the maximum and minimum average 
entanglement, respectively Ea and Ep. 

However, the measurement apparatuses that lead to 
some decompositions are not compatible with a continu- 
ous and local measurement process of the environments, 
or, as defined in f3|, the corresponding ensembles are not 



physically realizable through the continuous and local 
monitoring of the independent reservoirs. 

If one is to monitor the reservoirs continuously, then 
the evolution of the system respects a new equation where 
at each time step the system is projected onto a new 
state. Assuming that the system is initially prepared in 
a pure state, the evolution is then given by a sequence of 
pure states {|cro), \adt), |cr„dt)} which are called quan- 
tum trajectories of the system [1]. This process is also 
described by a CP-map in the form of ^ but now there 
is an extra restriction imposed over the different pos- 
sible unravellings K. Physically, this is related to the 
fact that some K's represent operations that combine 
measurements performed on the reservoir on different 
times, while the Markovian approximation that leads to 
the master equation imposes reservoirs with no memory 
which, under constant monitoring, must produce mea- 
surement outcomes that are not correlated in time. Fur- 
thermore, by imposing local observation of the reservoirs 
that act in each qubit, one also creates an extra restric- 
tion on the class of unitaries U that combine the different 
unravelling operations to unitaries of the type U^. ® l^B ■ 

Note that, in the limit of continuous monitoring there 
is yet another class of possible decompositions which is 
associated to the addition of classical noise to the mea- 
surement records (as in homodyne detection, for exam- 
ple [12 ). In this case, the decomposition reads G_|_^ — 
^ [SI ± UK.J] , where * indicates the absence of Kq and 

Go, and Go = 1 - i E,=i[GI+G,+ + gI_G,^]. The vec- 
tor $7 (whose elements are complex constants) represents 
the external entity that adds the extra noise. When |r2| 
tends to infinity we reach the quantum state diffusion 
limit ^ . When restricted to U = Uj^ ® Ub this decom- 
position always decreases the average entanglement as 
should be expected given that the added noise blurs the 
outcomes of the measurements on the reservoirs. This 
fact is used in to minimize the average concurrence. 

We show bellow that the local and continuous restric- 
tions made on the monitoring protocols of the reservoirs 
impose limitations on the minimum and maximum av- 
erage entanglement E attainable from realizable ensem- 
bles. Here, _£($) = '}2iPi-^i\'t'i))i where the entangle- 
ment of realizable state that happens with prob- 
ability Pi — {4>i\4>i), is obtained by calculating the Von 
Neumann entropy S{p) — — trjplogp} on the reduced 
density matrix over either one of the qubits, E{\(l)i)) — 

5(trB{|<^,)(0,|})- 

Let us first consider the example of both qubits ini- 
tially prepared in a maximally entangled state (for ex- 
ample one of |$±) — ^(|00) ± |11)) and subjected to 
independent dephasing reservoirs represented by Kraus 
operators given hy Aq — {\ — and Ai = Z^J^dt in 

the short time expansion (up to first order in dt). Here 
Z is the usual Pauli matrix, 7 is the system-environment 
coupling and the same holds for observer B. This partic- 
ular unravelling (Wa(b) = 1) already gives the E^ of the 



^B) 

system in each time interval dt 



Dephasing implies 
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that the system stochastically switches between |$+) and 
|$_) and by monitoring the reservoir one knows if the 
phase flip (the jump) has occurred or not. This result 
can still be classically communicated between the parts 
which then completely preserves the entanglement of the 
system. If, however, the observer is capable of phase flip- 
ping its local qubit, then not even the communication is 
required and the protocol becomes entirely local, similar 
to the one presented in [T3] for spontaneous emission. 
This also means that in this case EA[p{t)] corresponds to 
a locally realizable ensemble. Changing the unravelling 
also changes the decompositions and E, however the Ep 
of the unmonitored system is unreachable as shown in 
Fig [T] for the optimized local unravelling that minimizes 
E. 



still gives the same average values for all unravellings, as 
shown in |4l [15] . 
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FIG. 1: (Left) Dephasing reservoir for a Bell state: there is 
a gap between the average entanglement unraveled by local 
quantum state diffusion and the entanglement of formation 
Ef- The entanglement of assistance Ea and the maximum 
unraveled entanglement (MAXU) coincide. (Right) Dissipa- 
tive reservoir for a state presenting sudden death of entan- 
glement: quantum state diffusion still produces the minimum 
average entanglement which never goes to zero. In this case, 
there is a gap between both Ep and Ea, and the respective 
minimum and maximum locally realizable entanglement. Ea 
is in between the fidelity bound (F) and the eigenstate bound 
(E). 



Turning to dissipative environments, the U = TL un- 
ravelling features Kraus operators of the type Aq = 



1 - ^L^Ldt and Ai 



with L = |0)(1|. For ini- 



tial states in the subspace {|01), |10)} like, for example, 
= (|01) -I- |10))/-\/2, the entanglement of assistance 
of p(t) can once again be exactly obtained and coincides 
with the fidelity bound [7^ . This derives from the partic- 
ular time evolution of the system which leads to a mix- 
ture of a maximally entangled state and a separable one, 
p{t)=p{t)\^+){^+\ + [1 -p(t)]|00)(00|, and EA[p{t)] is 



given by p{t) 



-ft 



This maximum is exactly the one 



obtained by the Ua{b) = 1 unravelling [THl [H], which 
means that once again, local strategies allow for the max- 
imum extraction of entanglement from the system. In 
this case, classical communication is required given that 
a single decay completely destroys the entanglement of 
the system. Also note that, in this case, it is particularly 
important to use Ep as a lower bound for the average en- 
tanglement since the Wootters concurrence, for example. 



From the point of view of the detection of unreal- 
izable ensembles, the most interesting example, how- 
ever, appears when we analyze initial states that present 
finite time disentanglement such as \a) = ^(|00) -I- 

\/7|ll)) [16]. In this case, neither extremes of the en- 
tanglement of the system are achievable through local 
measurements on the independent reservoirs as shown in 
Fig. lb. For the Entanglement of assistance, this derives 
immediately from the fact that local strategies cannot 
increase the average entanglement while Ea is above the 
initial entanglement for times of the order of I/7. The 
most significant difference, however, is between the real- 
izable entanglement and Ep[p(t)] because the local mon- 
itoring of the reservoirs always preserve some entangle- 
ment in the system while for the unmonitored evolution 
entanglement disappears in a finite time |16j . 

The persisting gap between E and Ep necessarily 
means that ignoring the information that goes to the 
reservoir (master equation evolution of p) is different 
from continuously acquiring some of this information 
and then erasing it. The entropy of the measurement 
outcomes of a general unravelling is given by S{G) = 
-J^iPi^ogipi), with Pi = (GIGi) being the probabil- 
ity of the zth outcome. The information acquired from 
the reservoir is maximum if the entropy goes to zero, 
and minimum if the entropy is maximum. For exam- 
ple, for a single reservoir, the U — t unravelling in 
the short time expansion is given by an operator Ki 
proportional to ^/dt (jump operator) and its comple- 
mentary Kq = 1 — ^k\Ki. The entropy, given by 
5(K) = -pi log(pi) - (1 -pi)log(l -pi), with pi cx dt, 
tends to zero for small time steps — > 0, maximizing 
then the information. 

This unravelling locally produces the entanglement of 
assistance both for dephasing and spontaneous emission 
in the {|01), |10)} subspace. In both cases, the most prob- 
able event (no jump) is correlated to preserving entan- 
glement in the system. In fact, for dephasing, both click 
and no click produce equally entangled states, hence the 
Ea is locally realizable. Whereas, for the decay, no click 
still preserves entanglement, however, a click in the reser- 
voir kills entanglement completely, hence Ea decays. On 
the other hand, any other unravelling that combines the 
jump and no jump operators, necessarily increases the en- 
tropy, decreases the extracted information and mixes the 
above mentioned correlations. In the worst case scenario 

-1/^2' 

I/V2 1/^/2 

the information is lost and the average entanglement is 
minimized. For dephasing, it means that an event (either 
click or no click) will be equally correlated to both \^±) 
states, while for the {|01), |10)} decay, it will be equally 
correlated to the initial state and to |00). These limits 
establish the local unravelling bounds and any other lo- 
cal unravelling (more general Ws) will produce average 



of a 50 — 50 superposition {U = 



) all 
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entanglement in between them. 

A similar analysis can be carried out in the diffusive 
limit. The measurement operators can be described by 
Dj = ifo + G*J^ with G* = KM'^ and 3dt = {GJ + 
Gl)'^dt + dW, such that dWi are infinitesimal Wiener 
increments (gaussian distributed variables of variance dt) 
and dW'^dW = tdt. The short time expansion of the 
map that reproduces the density matrix is recovered with 

p{dt)^ Jd^i{3)DJaDl (4) 

with d/i(J) being the normalized gaussian measure over 
the measurement outcomes. In this case, the addition 
of classical nome,dWi, blurs the outcomes and destroys 
the acquisition of information from the reservoir, which 
leads the diffusive limit also to the minimum local av- 
erage entanglement, as we have shown in Fig. [l] of the 
paper. Once again, all the information is erased and the 
realizable entanglement reaches its minimum but not the 
EfW)]- 

Before concluding let us comment on the possible 
changes related to lifting the local observation restric- 
tion. That means allowing for operations that coher- 
ently superpose the results observed in each independent 
reservoir. In Fig. 1, we show that this strategy may also 
increase the average entanglement of the system. This 
happens to be the case for initial states presenting finite 
time disentanglement for spontaneous emission. In this 
case, both the no jump trajectory and the superposed 
jumps are correlated to increasing entanglement for the 
first part of the evolution (times of the order of I/7). 
While the no jump trajectory performs the optimum sin- 
glet conversion protocol [131 HZ], correlated jumps of the 
form G± = ^ [Ai (X) 1 ± 1 (g) Bi] produce Bell states for 
both qubits, increasing then, the average entanglement of 



the system. This strategy is still not enough to reach the 
Ea of the unmonitored evolution, but for certain times 
it redefines the lower bound for this quantity. In [T5] the 
authors use an equivalent strategy plus classical noise to 
reach the other extreme, continuously reproducing the 
concurrence of the unmonitored system. 

We have studied the entanglement that can be pro- 
duced in quantum open systems when the environments 
coupled to the system are continuously monitored. We 
have shown a gap between the locally realizable entan- 
glement and the extremes entanglement of assistance and 
formation of the unmonitored system. We have also re- 
lated the minimum attainable entanglement to the era- 
sure of the information provided by the readout of the en- 
vironment. On its turn, the above mentioned gap means 
that erasing this information is fundamentally different 
from not collecting it. This fact is better exemplified 
when the system is initially prepared in a pure entan- 
gled state and presents finite time disentanglement un- 
der dissipation, in which case it is impossible to locally 
unravel an ensemble with zero average entanglement at 
finite time. There is always at least one entangled tra- 
jectory which means that the average entanglement can 
only go to zero asymptotically in time. This fact in turn 
provides a method to protect entanglement and avoid 
entanglement sudden death. 
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